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Abstract. We investigate the slow-roll inflation in the nonminimal derivative coupling (NDC) 
model with exponential, quadric, and quartic potentials. It was known that this model 
provides an enhanced slow-roll inflation induced by gravitationally enhanced friction even for 
a steep exponential potential. In the phase portrait, the inflationary attractor is described 
by the slow-roll equation. Introducing the autonomous form, the inflation is regarded as an 
emergence from the saddle point and it leaves this fixed point along the slow-roll equation. 
We show explicitly that if one uses the NDC with potentials, the slow-roll inflation is easier 
to be implemented than the canonical coupling with the same potentials. 
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1 Introduction 

The nonminimal derivative coupling (NDC) [1, 2] was achieved by coupling the inflaton 
kinetic term to the Einstein tensor such that the friction is enhanced gravitationally at 
higher energies [3]. This gravitationally enhanced friction mechanism is a powerful tool to 
increase friction of an inflaton rolling down its own potential without introducing new degrees 
of freedom unlike hybrid inflation. Actually, this NDC makes a steep (non-flat) potential 
adequate for inflation without possessing higher-time derivative terms (ghost state) [4, 5]. 
It is well-known that the exponential potential provides a power-law inflation and thus, it 
cannot account for an inflationary theory since the inflation never ends and an additional 
mechanism is required to stop it [6] . The kinetic coupling flattens the potential effectively as 
well as it increases friction. Recently, this coupling allowed inflation to take place for a wide 
range of values A in the (non-flat) exponential potential V = Voe~^ which implies that it 
features a natural exit from inflation [7]. 

In this paper, we study the autonomous dynamical system of an homogeneous and 
isotropic configuration of a scalar field non minimally coupled to gravity, as in [8-16]. We 
found that in the two dimensional dynamical system of the field value and its derivative, 
the space of initial conditions providing inflationary attractors is larger in NDC than in the 
canonical case. Specifically we have investigated NDC models with exponential, quadric, and 
quartic potentials since non-flat potentials of exponential and quartic potentials with canon¬ 
ical coupling (CC) [17] are in tension with the Planck data [18]. Especially, the exponential 
potential is regarded as a testing potential because it could provide a power-law (eternal) 
inflation in the CC model while it could provide a slow-roll inflation in the NDC model [19]. 


2 An inflation model of NDC 

Let us consider an inflation model whose action includes NDC of scalar field </> with a potential 
term [3-5, 20] 



( 2 . 1 ) 
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where Mp is a reduced Planck mass, M is a mass parameter, and G fiv = R pu — g pu R/2 
is the Einstein tensor. Here, we do not include a canonical coupling (CC) term like as a 
conventional combination of (g pi , — G flll /M 2 )d^(j)d l '4 > [19, 21] because this combination won’t 
make the analysis transparent (see Appendix A for the autonomous system for CC+NDC). 

Varying the action (2.1) with respect to the metric tensor g leads to the Einstein 
equation 


G 


fllS - 


1 rpNDC 

M 2/J " 


( 2 . 2 ) 


where T^ DC takes the complicated form 


tiNDC 

1 fllS 


1 - 2V p 0V (p 0 J R£ ) + \G^(y<t>f 


M 2 L2 


-R 


jlpVCT 


V p </>V ff 4 


-V p V>V„V p </> + (V P V ^)V 2 0 


-g^[ - R pc y P <t>v^ + ^(v 2 </>) 2 - i(v"vv)v P v c 


(2.3) 


Here, we note that even though fourth-order derivative terms are present in T^ DC , there 
is no ghost state which means that any higher-time derivative term more than two is not 
generated 1 . On the other hand, the scalar equation is derived to be 

--t-G^V^^-V' = 0, (2.4) 

M z 

where the prime (') denotes derivative with respect to qb. 

In this work, we consider a spatially flat spacetime by introducing cosmic time t as 


ds 2 = g pi ,dx p dx v = — dt 2 + a 2 {t)8ijdx l dx\ 


(2.5) 


where a(t) is a scale factor. In this spacetime, two Friedmann equations and scalar equation 
become 


H 2 = 


3 M 2 
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( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


where H = a/a is the Hubble parameter and the overdot Q denotes derivative with respect 
to time t. We observe from (2.6) that the energy density for the NDC is positive (ghost-free) 
[ 22 ]- 

Now we define two slow-roll parameters as 


CN = 


H 
~H 2 ’ 


S K = A 

H(j) 


(2.9) 


1 We check easily that the two terms for <f> 2 in the second line of Eq.(2.3) cancel against each other and the 

last two terms in the last line of Eq.(2.3) do too. 
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Imposing slow-roll conditions (e^ < 1, 5j \r <C 1), Eqs. (2.6)-(2.8) can be written approxi¬ 
mately as 


H 2 ~ 


3 M 2 


V, 


H ~ - 


3 H 2 

2MpM 2 ' 
M 2 

3H * “ * 3iP V ' 


j.2 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


We note that the slow-roll parameter ejy can be written as 


V H 

e N ~ e N ~ e N , 


where the potential ejy and Hubble slow-roll parameters ejy are given by [4, 23] 


v MpM 2 V' 2 

e N = 


6 V 2 H 2 ’ 


e N = 


2 MpM 2 


(2.13) 


(2.14) 


The end of inflation can be identified with the end of slow-roll regime. Therefore, we may 
define the value of field at the end of inflation 0/ by making use of = 1. 


3 Dynamical analysis 

In this section, we wish to perform the dynamical analysis for (2.6)-(2.8). For this purpose, 
one introduces with some dimensionless quantities and obtains their autonomous system. 
Before we proceed, it is instructive to investigate the canonical coupling (CC) model for 
comparison and exercise. 


3.1 Analysis of CC models 

In this case, the corresponding field equations are given by 


H 2 = 


3 M 2 


1 V 


T V 


H = - 


2 M 2 


4> + 3H<j> + V' = 0. 


(3.1) 

(3.2) 

(3.3) 


It is well-known that in the slow-roll approximation, they reduce to 


9 V 
H 2 ~ 


3Mp ’ 


H = — 


2 M 2 


3H(j) ~ -V', 


(3.4) 

(3.5) 

(3.6) 


which are obtained by imposing slow-roll condition of ec = — H/H 2 1 and 5c = 4>/H4> <S 

1. Here we have ec — — e c 


V — Mp V 2 h _ 

e ° = — V 2 ’ €c = 2MpH 2 ' 


(3.7) 
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We note that the NDC parameters of e^ H (2.14) are slightly different from the CC parameters 
e^n (3.7). In particular, e]v = e c x M 2 /(3H 2 ) implies that when M 2 /(3H 2 ) <C 1, inflation 
in the NDC can happen more easily than the CC case. 

To perform the dynamical analysis for the cosmological evolution equations (3.1)-(3.3), 
we consider the dimensionless variables 


x = 


MpV$H' 


W _ ^ V ' 

-, ay = v6Mp —, 

MpVZH V 


r = 


VV" 

-yTT- 


(3.8) 


The connections between ( x,ay) and (e^,e^) are given by 

x 2 = € f, a 2 v = 12e^. (3.9) 

Then, the first Friedmann equation (3.1) becomes a constraint equation 


2 , 2 

x + y z 


= 1 , 


(3.10) 


which implies that the variable y can be eliminated from the dynamical equations. 

It is found that (3.1)-(3.3) lead to the autonomous form for X = ( x,ay ) 

3(x 2 — l)(x + (3.11) 

o^x(r-l) (3.12) 

with N = In a. Fixed points of first-order autonomous system are given by dx/dN = 0 and 
day/dN = 0. By a fixed point we mean that x and ay do not change as the universe evolves. 
Perturbing these points leads to stable point (attractor), saddle point, and unstable point 
(repeller). Then, one can classify these according to their eigenvalue solution to perturbed 
equations: all negative (stable point); all positive (unstable), different signs (saddle point). 

At this stage, it is worth noting that the slow-roll trajectory can be found by considering 
two different ways: One is found by solving the slow-roll equations (3.4)-(3.6) and the other 
is found by solving the autonomous system (3.11) and (3.12) in the slow-roll approximation. 
The former indicates an inflationary attractor in phase portrait (</>, </>), while the latter shows 
that inflation is regarded as “an emergence from the saddle point and it can leave this 
fixed point along the slow-roll equation”. For this purpose, we introduce three potentials: 
exponential, quadric (chaotic), and quartic potentials. The first one is regarded as a testing 
potential because it cannot provide a slow-roll inflation in the CC model, whereas it could 
provide a slow-roll inflation in the NDC model. 

We should distinguish between inflationary attractor in (0, 0) and attractor of stable 
fixed point in (x,ay). To this end, we note that slow-roll equations (3.4)-(3.6) are combined 
to give 

<P-~Mp^-, (3.13) 

which describes an inflationary attractor. This can be expressed in terms of x and ay defined 
in (3.8) as 

x ~ — —. (3-14) 


dx 

dN 

day 

~dN 
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Figure 1. The phase portrait (cf>,(j)) [left] and stream flow (x,y) [right] for V = Voe A< ^ with Vo = 0.1, 
A = 1 and Mp = 1. Left panel shows the trajectories (green curves) with various initial conditions 
toward the inflationary attractor (red curve). This red curve is mapped into a fixed point of red dot 
at (—0.408,0.912) in the right panel. This corresponds to a stable fixed point (attractor), while the 
point (0,1) denotes a saddle point. The green curves in the left are mapped on the green curves on 
the upper circle in the right. 


We note that Eq.(3.14) satisfies e// ~ which yields x 2 — ay/36, when substituting (3.9) 
into . Also, (3.14) describes a slow-roll line solution to an approximate equation 


dx 

dN 


, ay n 
~ X + — ~ 0 , 

fa 


(3.15) 


which was found from (3.11) by taking into account the slow-roll condition of ef/, = 3x 2 <C 1. 
This implies that any fixed point can be defined only for T = 1 (ap=const) or x = 0 (ay=0), 
while the slow-roll line is defined without imposing cq/ =cons t- Here, a slow-roll line is given 
by a function [ay{x) = — 6x\ which exists during inflation. Hence, we identify the inflationary 
attractor in phase portrait (</>,</>) with the slow-roll line (equation) in stream flow (j,«y). 
The slow-roll approximation reduces the order of the system equations by one and thus, its 
general solution contains one less initial condition. It works only because the solution to the 
full equations has an attractor property eliminating the dependence on the extra parameter. 


For complete analysis, we choose an explicit potential, 
(z) V = H o e A0 


For an exponential potential, the inflaton velocity (3.13) is given by 


<j> — —MpA 



(3.16) 


Solving (3.1)-(3.3) numerically for various initial conditions, they show the trajectories. Also, 
the slow-roll equation (3.16) may indicate an inflationary trajectory [see Fig.l (left)]. How¬ 
ever, this case corresponds to = Mp/2 (power-law inflation), implying that it cannot be 
a complete inflationary model because inflation never ends with the potential V = . 

In order to see it more clearly, we have to know what happens in the stream plot where 
fixed points are included naturally. In this case, we have T = l(ap=const). Therefore, we 
could not make a stream flow on (i,ay). Instead, we have the stream flow on (x,y). We 
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Figure 2. The phase portrait (</>, <p) [left] and stream flow (x,ay) [right] for V = Vo</> 2 with Vo = 0.1 
and Mp = 1. Left panel shows the trajectories toward the inflationary attractor (red lines) for cf> ^ 0 
including a spiral sink. These red lines correspond to a red line of ay = —6x in the right panel. Also, 
the green flows in the right show the numerical plot for (x,ay) when using (3.1)-(3.3). These flows 
indicate that the inflation is emergent from the saddle point (•) (0,0) and is realized along the red 
line. Finally, the inflation ends at the point (■) (=f 0.46, ±3.48), corresponding to = 1. 


find a corresponding fixed point of ( x , y ) = (—ay/6, yj 1 — ay/ 36) = (—0.408, 0.912) in Fig.l 
(right) which turns out to be an attractor (stable fixed point), in addition to the (0,l)-saddle 
point. Clearly, it indicates that the potential V = Vfoe ^ in the CC model cannot provide a 
complete slow-roll inflation because power-law inflation never ends. This means that there 
is no spiral sink which indicates the end of inflation in the left of Fig. 1. 

(ii) V = V 0 </> 2 


In the case of chaotic potential, the inflationary attractor (3.13) is given by 


f-2M PV /f, (0> 0) 

\+ 2 MpJf, {$< 0 ) 


(3.17) 


which corresponds to a positive (negative) constant (/>-line for 4> < 0 (0 > 0) in the (0, <j>) 
picture. 

Figure 2 shows a typical picture based on numerical computation. When the universe 
evolves according to (3.1)-(3.3) and (3.11)-(3.12), the inflationary attractor and slow-roll line 
are given by (3.17) in the phase portrait (0, <fi) [Fig.2 (left) panel] [6, 13, 14] and (3.14) in 
stream flow ( x,ay ) [Fig.2 (right) panel], respectively. There are three phases in the CC 
case [24]: initially, kinetic energy dominates and due to the rapid decrease of the kinetic 
energy the trajectory runs quickly to the inflationary attractor line (3.17). All initial trajec¬ 
tories are attracted to this line, which is the key feature of slow-roll inflation. Finally, at the 
end of inflation, there is inflaton decay and reheating (spiral sink). 

On the other hand, the inflation can be realized as an emergence from the potential- 
dominated fixed point (0,0) which is a saddle point because orbits near it are attracted along 
one direction and repelled along another direction. The stream flows in the right of Fig. 2 
can leave this fixed point only along the red line of ay = — 6x corresponding to the slow-roll 
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Figure 3. The phase portrait (<f>, <fi) [left] and stream flow (x, ay) [right] for V = Vo^> 4 with Vq = 0.1, 
Mp = 1. Left panel shows the trajectories toward the inflationary attractor (red lines) followed by 
spiral sink. These lines correspond to a red line of ay = — 6x in the right panel. Also, the green flows 
in the right show the numerical plot for ( x , ay) when using (3.1)-(3.3). These flows indicate that the 
inflation is emergent from the saddle point (•) (0,0) and is realized along the red line. Finally, the 
inflation ends at the point (■) (±0.50, ±3.48), which corresponds to = 1. 

line (3.14) and the inflation ends at the point (±0.46, ±3.48), which corresponds to = 1 
with ay = ±3.48. We note here that cc-axis is the stable manifold of the saddle point, while 
ay = —6x is the unstable manifold of the saddle point. 

(in) V = V 0 c/> 4 

The inflationary attractor (3.13) for this case is given by 


0 





(3.18) 


Figure 3 shows that when solving (3.1)-(3.3) and (3.11)-(3.12) numerically, the inflationary 
attractor (decreasing function) is given by the slow-roll equation (3.18) followed by a spiral 
sink in the phase portraits (0, </>) [Fig-3 (left)] and the red line (3.14) in stream flows ( x,ay ) 
[Fig.3 (right)], respectively. Inflation can be realized as an emergence from the saddle point 
(0,0). The inflation can leave this fixed point only along the red line of ay = — 6x correspond¬ 
ing to the slow-roll line (3.14) and it ends at the point (±0.50, ±3.48), which corresponds to 
€q = 1 with ay = ±3.48. 


3.2 Analysis of NDC models 

Now we turn to the NDC case. To perform the dynamical analysis for (2.6)-(2.8), we consider 
the dimensionless parameters u and a u in addition to y and T defined in (3.8) 


13 

u = \t -- 


2 M P M 


a, 


= V6 


My MV' 
VH ' 


The connections between ( u,a u ) and (e^,e)y) are given by 

2 _ H 2 _ o/? V 

u — e N , a u — 36e N . 


(3.19) 


(3.20) 
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Then, the constraint equation (2.6) becomes 


v 2 + y 2 = 1 , 


(3.21) 


which implies that one may eliminate y from the dynamical equations. It turns out that 
(2.6)-(2.8) lead to the autonomous form for X = (u,a u ) 


du . . 1 . 9 , 

— = (2£at - 3 )u - -a u ( 1 -u), 


da. 


1 


d N = €NOlu + ^ ua «’ 


where e^v = —H/H 2 satisfies the relation 


1 


(1 + u 2 )eN = 3 u 2 + -a u u( 1 - u 2 ). 

O 


(3.22) 

(3.23) 

(3.24) 


On the other hand, we note that slow-roll equations (2.10)-(2.12) give us the inflaton velocity 

V' 


0 ~ — MpM 2 


3HV’ 


(3.25) 


which corresponds to 


u^-^. (3.26) 

b 

We note that Eq.(3.26) satisfies ~ e]y which yields u 2 ~ a 2 /36, when substituting (3.20) 
into e^’ . Also, it is worth to mention that Eq.(3.26) corresponds to Eq.(3.14) in CC. Fur¬ 
thermore, Eq.(3.26) can be realized as a slow-roll line solution to (3.22) when implementing 
slow-roll condition of = « 2 C 1 which takes the form 


du 

dN 


a„ 

~ — 3u -~ 0. 

2 


(3.27) 


Differing with the CC model, there exists an upper limit of cj) 2 


o < 4> 2 < 4>l = -MpM 2 , 


(3.28) 


which comes from Eq.(2.6) yielding H 2 { 1 — ^> 2 /^ 2 ) = V/3Mp where the left-handed side 
should be positive for a positive potential. This may be regarded as another representation 
of slow-roll condition (ejy = ii 2 < 1 ), 

Now, we solve (2.6)-(2.8) and (3.22)-(3.23) numerically by considering three of expo¬ 
nential, quadric (chaotic), and quartic potentials. 


(*) V = VoeM 


For an exponential potential, the inflaton velocity (3.25) is given by 






M^M 2 A _a 
— , e 2 

V3W 


(3.29) 
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Figure 4. The phase portrait {<f>, <p) [left] and stream flow (u, a u ) [right] for V = Voe Xcl> with Vq = 0.1, 
A = 1, Mp = 1, M = 0.1, and (f) c = MpMy/2/\f3 ~ 0.08. Left panel shows the trajectories toward 
the inflationary attractor (red curve). Here is no stable limited cycle. This curve corresponds to a 
red line of a u = —6 u in the right panel. Also, the green flows in the right show the numerical plot for 
(■ u,a u ) when using (2.6)-(2.8). These flows show that inflation is emergent from the saddle point (•) 
(0,0) and is realized along the red line. Finally, the inflation ends at the point (■) (—0.69,6), which 
corresponds to = 1. 


which implies that when solving (2.6)-(2.8) numerically for various initial conditions, any 
solution should attract the trajectory of (3.29) [see Fig.4 (left)]. Here we note that there is 
no spiral sink because the potential approaches zero in the limit of <f> — oo. We point out 

that a u = —6 u (3.26) cannot be an attractor (point) because of evolution of H in a u (3.19) 
even though V'/V = A. This contrasts to the CC case, which yields an attractor (point) at 
x = —av/Q for the exponential potential [see Fig.l (right)] implying that the slow-roll never 
ends. 

On the other hand, for the NDC case with the exponential potential, one checks that 
the inflationary attractor in the ( u,a u ) is given by a slow-roll line of a u = —6 u [see Fig.4 
(right)]. In this case, inflation can be realized as an emergence from the saddle point (0,0). 
The stream flows in the right of Fig. 4 can leave the point (0,0) only along the red line 
of u = —a u /6 corresponding to the slow-roll line (3.26) and the inflation ends at the point 
(—0.69, 6), which corresponds to = 1 with a u = 6. 

(ii) V = Vq cf ) 2 


In case of chaotic potential, the inflaton velocity (3.25) is given by 


<j> ~ 


2AfjjM 2 i 
~ x/3Vo <F’ 
, 2 M^M 2 i 

+ vm 0 *’ 


(0>O) 

(0 < o) 


(3.30) 


Figure 5 shows that when solving Eqs. (2.6)-(2.8) and (3.22)-(3.23) numerically, the infla¬ 
tionary attractors are given by the red curve (3.30) followed by stable limited cycle on (0, </>) 
[Fig.5 (left)] and the slow-roll line (3.26) starting at (0,0) [Fig.5 (right)]. We note here that 
the stable limited cycle where nearby curves spiral towards closed curve C appears instead 
of the spinal sink in the CC. Hence, the trajectories trace out a closed curve C. 

Inflation can be realized as an emergence from the saddle point (0,0). The stream flows 
in the right of Fig. 5 indicate that the inflation can leave the saddle point (0, 0) only along 
the red line of u = —a u /6 corresponding to the slow-roll line (3.26). The inflation ends at the 
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Figure 5. The phase portrait (cj>, </>) [left] and stream flow (u, a u ) [right] for V = Vq(/) 2 with Vo = 0.1, 
Mp = 1, M = 0.1, and <f> c = MpMy/2/\/3 ~ 0.08. Left panel shows the trajectories toward the 
inflationary attractor (red curves) followed by stable limited cycle. These curves correspond to a 
red line of a u = —6 u in the right panel. Also, the green flows in the right show the numerical plot 
for ( u,a u ) when solving (2.6)-(2.8). These flows indicate that the inflation is emergent from the 
saddle point (•) (0,0) and is realized along the red line. Finally, the inflation ends at the point (■) 
(yO.63, ±6), which corresponds to ej) = 1. 


point (t 0.63,±6), which corresponds to (X; = 1 with a u = ±6. Furthermore, comparing the 
red curve in Fig.5 (left) with the red line in Fig.2 (left) indicates that the NDC suppresses 
the kinetic term tf> whereas it enhances higher friction. 

(in) V = V 0 cf > 4 


The inflaton velocity (3.25) for this case is given by 


4MpM 2 1 


(3.31) 


Finally, Figure 6 shows that when solving (2.6)-(2.8) and (3.22)-(3.23) numerically, the infla¬ 
tionary attractors are given by the red curve (3.31) followed by stable limited cycle on (0, </>) 
[Fig.6 (left)] and slow-roll line (3.26) on ( u,a u ) [Fig.6 (right)]. Inflation can be realized as 
an emergence from the saddle point (0,0). The stream flows in the right of Fig. 6 indicate 
that inflation can leave the point (0,0) only along the red line of a u = —6 u corresponding to 
the slow-roll line (3.26) and the inflation ends at the point (=p0.66,±6), which corresponds 
to fXj = 1 with a u = ±6. Also, comparing the red curve in Fig.6 (left) with the red curve 
in Fig.3 (left) indicates that the NDC suppresses the inflaton velocity (j) whereas it enhances 
higher friction. This means that the initial kinetic energy in the NDC is small when com¬ 
pared with the potential energy. On the other hand, the kinetic energy in the CC rapidly 
decreases until eventually the potential energy dominates and the universe may enter the 
slow-roll inflation. However, it is not easy to obtain an enough e-folds number from 0 4 in 
the CC. This explains why this potential was ruled out [18]. 


4 Summary and Discussions 

First of all, we compare the NDC with the CC with the same potentials. In the phase portrait 
(4>, cf >), the slow-roll inflation is identified with the presence of inflationary attractor (slow-roll 
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Figure 6. The phase portrait (cj>, </>) [left] and stream flow (u, a u ) [right] for V = Vq(/) 4 with Vq = 0.1, 
Mp = 1, M = 0.1, and cf> c = MpMy/2/y/3 ~ 0.08. Left panel shows the trajectories toward the 
inflationary attractor (red curves) followed by stable limited cycle. These curves correspond to a red 
line of a u = —6 u in the right panel. Also, the green flows in the right show the numerical plot for 
( u,a u ) when using (2.6)-(2.8). This flow shows that the inflation is emergent from the saddle point 
(•) (0,0) and is realized along the red line. Finally, the inflation ends at the point (■) (=f 0.66,±6), 
which corresponds to ej) = 1. 


equation) followed by spiral sink (0,0) (the ending of inflation) for the CC and stable limited 
cycle centered at (0,0) for the NDC. See Fig. 7 for the different behaviors of 4> and 4> between 
the CC and NDC during the evolution of the universe. On the other hand, the corresponding 
autonomous forms for X=(x(a),ay( u )) indicated that the slow-roll inflation is realized by 
emerging a saddle point (0,0) (the beginning of inflation) and inflation leaves this point along 

the slow-roll equation. This is so because of ( x,ay ) — (ln[<^>], ^) for cj) 2 and (^, ^) for <p A in 
the CC case. In the NDC case, ( u,a u ) ~ (cj), ^) and (<p, ). This implies that even though 

the limit of <f> —>• 0 is not implemented by the autonomous system, the other limit of <f> —>• oo 
can be easily realized at the origin (0,0). 

We summarize our results in Table 1. 0* is the value of cj) when the phase space trajectory 
hits the inflationary attractor while 0/ is determined by N = 1. The relations are given 
for the exponential (e^) and power-law potential ( <j ) n ) by 

= —XMpN (CC), e A< ^ - e A ^ = -A 2 M^M 2 ^ (NDC) 

Vo 

$ -$ = —2MpnN (CC), <^+ 2 - £ +2 = -n{n + 2)M^M 2 ^ (NDC), 

j Vo 

where the e-fold number is IV = 70. Acfi = cj)f — fa in the CC and NDC are 15.38 and 2.07 
for cj) 2 , 21.01 and 1.39 for </> 4 , and not available (N.A.) and 4.94 for e A ^. During inflation, the 
inflaton must be in a slow-roll regime. This implies that the inflaton varies little during the 
inflationary phase and thus, it satisfies f <7 and <j) ^ 3 H<j). In this sense, the NDC with 
potentials is easier to make slow-roll inflation than the CC with the same potentials. 

Especially, the exponential (non-flat) potential is regarded as a testing potential because 
it provides a power-law inflation in the CC model while it could provide a slow-roll inflation in 
the NDC. We note that as was shown in Fig.l (right), there is no saddle point as a starting 
point of slow-roll line even though (0,1) denotes a saddle point. Also, there is no stable 
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Figure 7. The evolution of <j>(t) [left] and <j>(t) [right] with respect to time t for potential <jr. The 
left figure shows that the inflaton varies little during large inflationary period for the NDC, while it 
varies quickly during small inflationary period for the CC. After inflation, <j> decays with oscillation for 
CC, while it oscillates rapidly for NDC. The right one indicates that for large t, <j> oscillates without 
damping for the NDC, while it oscillates with damping for the CC. 



CC 

potential 

(pi 

<t>f 

spiral sink 

saddle point 


16.79 

1.41 

0 

0 


23.83 

2.82 

0 

0 

e x<t> 

N.A. 

N.A. 

X 

X 


NDC 

potential 

(pi 

<t>f 

stable limited cycle 

saddle point 

4> z 

2.74 

0.67 

0 

0 


2.35 

0.96 

0 

0 


1.95 

-2.99 

X 

0 


Table 1. Summary on how inflation takes place for the CC and NDC with three potentials of </> 2 , <^> 4 , 
and e A< A <pi denotes minimal value satisfying the e-folds N > 70 and <f>f is the end of inflation. 
The presence of either spinal sink or stable limited cycle represents the end of inflation in the phase 
portrait (</>, </>), while the saddle point denotes the emergence of inflation in the autonomous form. 


limited cycle in the phase portrait [Fig.4 (left)] of the NDC which features the exponential 
potential. 

In order to see how slow-roll inflation occurs in the NDC case, we compare the NDC 
with the CC when one chooses the chaotic potential V = Vq< p 2 . Figure 7 (left) indicates that 
the inflaton varies little during large slow-roll period for the NDC, while it varies quickly 
during small slow-roll period for the CC. This implies that the NDC would inflate enough 
comparing with the CC. After the inflation [see Fig.7 (right)], the velocity of the inflation (f> 
oscillates with respect time in the NDC, while its velocity describes a damped-oscillation in 
the CC. This distinguishes the stable limited cycle for the NDC from the spinal sink for the 
CC. 

It is worth noting that there are three phases (—> • —^-spiral sink) in the CC case [24]: 
i) Initially, kinetic energy dominates [see Fig.2 (left) and Fig.3 (left)], ii) Due to the rapid 
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decrease of the kinetic energy, the trajectory runs to the inflationary attractor line (3.17). 
All initial trajectories are attracted to this line, which is the key feature of slow-roll inflation, 
iii) At the end of inflation, the magnitude of inflaton velocity \(j>\ decreases. There is inflaton 
decay and reheating which shows the appearance of spiral sink. On the other hand, three 
stages (—> • —^-stable limited cycle) in the NDC are as follows: i) Initially, potential energy 
dominates [see Fig.5 (left) and Fig.6 (left)], ii) Due to the gravitationally enhanced friction, 
all initial trajectories are attracted quickly to the inflationary attractor, iii) At the end of 
inflation, the magnitude of the inflation velocity \4>\ increases. Then, the stable limited cycle 
appears, which differs from the spiral sink in the CC case. This is a rapidly oscillating phase 
[25], but it appears after slow-roll inflation [26]. 

We would like here to warn the reader that, because of numerical complexity, all our 
analysis are either made in the NDC or in the CC system but not in the more physical 
combined NDC+CC system (although some hints are given in the Appendix). Nevertheless, 
we expect the sector we have excluded to provide an even larger space of initial conditions 
leading to inflation. Thus, our results although limited already brings the important message 
that the NDC coupling facilitate inflation with respect to the sole CC one. 

Finally, we wish to mention that for the power-law potential <^> n , the field excursion of 
the inflaton is sub-Planckian due to the NDC [17]. Thus, the tensor-to-scalar ratio r is a 
factor of (n + 2)/2 smaller than the results in the CC which brings the quartic and quadric 
potentials to be consistent with the observation at the 95% CL. 
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Appendix 

A The autonomous system for CC+NDC 

In this appendix, we consider an action including a canonical coupling term as 

1 


S= drxyj—g 


4_ 


L 2 


9fiU m 2 Gflu 




(A.l) 


For a spatially flat spacetime (2.5) with 4> = (j>(t), the Friedmann equations and scalar 
equation are given by 


H z = 


1 


3 M 2 L2 


1 V 


1 + 


9 H 2 
M2 


+ K 


H = — 


2Mp 


1 + 


3 H 2 H \ 2 H , ■: 


H \ 2 H 

M 2 ~ M 2 ) ~ M2 ' 


1 + 


3 H 2 
M2 


3 H 2 2 H 


4> + 3 H (1 -|————b 




M 2 M 2 / 


+ V' = 0. 


(A.2) 

(A.3) 

(A.4) 
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To perform the dynamical analysis for the equations (A.2)-(A.4), we first introduce the 
following dimensionless quantities, 


x = 


MpV&H' 


u = 


2 M P M 


y = 


w 


M P V3H } 


ay 


= \/6Mp 


r 

v : 


r MpMV' k , 

= v/6 VH ■ (A.5) 


Here, x and u represent the dominance of the kinetic term for CC and the kinetic term for 
NDC, while y denotes the dominance of the potential term for CC and NDC. However, a 
constraint equation obtained from Eq. (A.2) is given by 


2,2,2 i 

x +u +y = 1, 


(A.6) 


which allows us to eliminate y. 

Introducing N = In a, (A.2)-(A.4) can be written in terms of the quantities (A.5) being 
consisting the following autonomous form for X = ( x,u,ay,a u ): 


dx . „ 

IN= X(C + S) ' 

du . 

Jn = uS ' 


day 

~dN 

da u 


= ayxiT - 1 ), 


1 


dN = eoiu+ 


(A.7) 

(A.8) 

(A.9) 

(A.10) 


where e, <5, and their relation are given by 


H A ^ 

e ~ ° i e = 

di 2 (j)H 


3x 2 + u 2 — | u 2 d 


(A.ll) 


It is not an easy task to analyze the autonomous system for X = (x,u,ay,a u ) (CC+NDC) 
completely because of its complexity. However, we expect to extract some information on 
the full system (CC+NDC) by analyzing fixed points obtained from the autonomous system 
(A.7)-(A.10) whose compact form is rewritten by X 7 = f(X). The fixed points X^ are 
extracted from the condition of X 7 = 0 and they provide qualitative information on the 
global dynamics of the system, independently of the initial conditions and specific evolution 
of the system. This information might include all previous fixed points and new fixed points 
from other combinations. First of all, we have recovered all previous fixed points from X 7 = 0 
and summarize them (Pl± ~P8) in Table 2. Then, what are new fixed points which might 
be found from other combinations? We expect to have three candidates from the analysis of 
X 7 = 0. It is checked easily that the two kinetic-dominant fixed point of (a;7,0,tt7,0,y7 ~ q) j s 
not allowed for the autonomous system (A.7)-(A.10) because both CC and NDC are kinetic 
term. The two remaining points are the potential-dominant fixed points of (0,ay,0,al,yf) 
and (0 ,ay,uf,al,yf). Explicitly, the former includes two saddle points P9 and P10, while 
the latter is given by repeller Pll± which is similar to P4-j-. Finally, it is worth to note that 
P9 and P10 might represent slow-roll inflation in the full system (CC+NDC) because these 
belong to saddle points. 
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fixed points 

dy , U , dw) £/) 

potential 

type 

stability 

where 

Pl± 

(±1,77 0 , 0 , 0 ) 


CC 

R 

Fig. 1 (right) 

P2 

(-1/7776,0,0,7^6) 

V 0 e* 

cc 

A 

Fig. 1 (right) 

P3 

(0,V6,0,0,1) 

V 0 e* 

CC 

S 

Fig. 1 (right) 

P4± 

(0,76,±1,0,0) 

V 0 e* 

NDC 

R 

Fig.4(right) 

P5 

(0,76,0,0,1) 

Voe* 

NDC 

S 

Fig.4(right) 

P6 

(0,0,0,0,1) 

v 0 r 

CC 

S 

Fig.2,3(right) 

P7± 

(0,0,±1,0,0) 

v 0 r 

NDC 

R 

Fig.5,6(right) 

P8 

(0,0,0,0,1) 

v 0 r 

NDC 

S 

Fig. 5,6 (right) 

P9 

(0,76,0,a u ,l) 

V 0 e* 

CC+NDC 

S 

new 

P10 

(0,G(\/ , 051 ^ 

v 0 r 

CC+NDC 

s 

new 

PH± 

(0,ay ,±1,0,0) 

v 0 r 

CC+NDC 

R 

new 


Table 2. List of all previous and new fixed points X^ in the autonomous system of CC+NDC 
for V = Voe^ and V = Vo<jP(jp = 2,4). In the stability, A, S, and R represent attractor, saddle 
point, and repeller, respectively where R means kinetic-dominant fixed point (y = 0) and S implies 
potential-dominant fixed point (y = 1). 
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